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We study the effect from a sequential fourth generation quark on penguin-dominated two-body 
nonleptonic B meson decays in the next-to-leading order perturbative QCD formalism. With an 
enhancement of the color-suppressed tree amplitude and possibility of a new CP phase in the elec- 
troweak penguin, we can account better for A C p{B° -> K + tv~) - A C p{B + -> K + tv°). Taking 
|Vt' a Vt'b| ~ 0.02 with phase just below 90°, which are consistent with the b — * s£ + £~ rate and the 
B s mixing parameter Arm, , we find a downward shift in the mixing-induced CP asymmetries of 
B° — > Kstt and <pKs. The predicted behavior for B° — » p°Ks is opposite. 

PACS numbers: 12.60.-i, 13.25.Hw, 12.38.Bx 
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CP violation (CPV) in b — ► s transitions is at the 
forefront of our quest to understand flavor and the 
origins of CPV, offering one of the best probes for 
New Physics (NP) beyond the Standard Model (SM). 
Several hints for NP have emerged in the past few 
years. For example, a large difference is seen in direct 
CP asymmetries in B — > Kir decays 

Ak-k = ^4cp(-B° 
A Kn o =A C p(B+ 



K 



= -0.093 ±0.015, 

= +0.047 ±0.026, (1) 



or AAkx = A Kn o - Ak-k = (14 ± 3)% 0- As it was 
not predicted when first measured in 2004, it has stim- 
ulated discussion on the potential mechanisms that 
may have been missed in the SM calculations 0, 0, [H ■ 
Better known is the mixing-induced CP asymmetry 
Sf measured in a multitude of CP eigenstates /. For 
penguin-dominated b — > sqq modes, within SM, S sq q 
should be close to that extracted from b — > ccs modes. 
The latter is now measured rather precisely, S cEs = 
sin 20i = 0.674±0.026 [y], where <pi is the weak phase 
in Vtd- However, for the past few years, data seem to 
indicate, at 2.6 a significance, 



AiS = S sq q 



Sees < 0, 



(2) 



which has stimulated even more discussions. 

Since the two modes in Eq. fl} differ by the sub- 
leading color-suppressed tree amplitude C and elec- 
troweak penguin amplitude -P EW , it is natural to con- 
jecture that C , -F^w, or maybe both have been un- 
derestimated (see [3] and references therein). In this 
paper we study the fourth generation effect on the 
B — > Kir, <j)K and pK decays in the perturbative 
QCD (PQCD) approach, which exploits both C" and 
P EW to account for Eq. ([I}, and gives the right ten- 
dency towards explaining Eq. (JSJ. 

At leading order (LO), PQCD predicted [3 the sign 
and strength of Ak-k-, but Ak-k was a l so negative and 



not too different from Akk ■ Going to next-to- leading- 
order (NLO) P, it was shown that inclusion of the ver- 
tex corrections increases C by a factor of 3, without 
affecting much the branching ratios, Akk becomes 
much closer to zero, although still negative. However, 
like many other studies in the literature 0, [l(| [HI], 
NLO PQCD also predicts AS Ks n° = S Ks k° —S C cs > 
within SM, opposite to what is indicated by Eq. @. 
Thus, if AS stands the scrutiny of time, one would 
still need to go beyond SM. 

The other path to account for Eq. {TJ), to have 
an effect coming from -P EW , in fact requires NP 
CPV. Among the available models, the sequential 
fourth generation is found to modify only P EW sig- 
nificantly [12j |. but not other amplitudes such as the 
QCD penguin P' . With reasonable parameters, e.g. 
with m t > ~ 300 GeV, the product of quark mixing 
matrix elements | V*' s Vt'b | = 0.01 ~ 0.03 with phase 
close to 90°, it was found in PQCD at LO that Akk 
could also be rendered vanishing. Remarkably, the 
constraints of the b — > s£ + £~ rate and the B s mi xing 
parameter Ams s are also satisfied. It was shown [1 31 ] 
in QCD factorization (QCDF) at NLO (and PQCD 
at LO for B — * Kir), that the same parameter choice 
produces a downward shift for Sk s k° an d S^,k s - The 
prediction is far from reliable, however, due to uncon- 
trollable hadronic uncertainties in QCDF. 

It is worthwhile, then, to reanalyze the fourth gen- 
eration effect using the NLO PQCD formalism, in part 
to find whether the preferred parameter space [12| is 
affected, and also to investigate the efficacy of combin- 
ing the C and -P EW approaches on CPV in penguin- 
dominated b — > sqq modes. It should be emphasized 
that both the strength and phase of V t 1 s Vfb would 
be effectively pinned down by the recent precise mea- 
surement of Amfl l [3], which will then lead to rather 
definite predictions for CPV in B s mixing, sin2$s s , 
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as well as D° mixing [l5|. We show that, using the 
same fourth generation parameters as before, we can 
explain the trend of the observed AS in B — » Kir 
and 4>K decays simultaneously, while understanding 
of AAk-k is also improved. We point out that the 
fourth generation effect is opposite in B — > pK, which 
increases S p ok s from a low value in SM [lOL Hq. 

It is instructive to elucidate the underpinnings of 
the two proposals to resolve the B — ► Kir puzzle, i.e. 
enhanced C by vertex corrections in SM, vs fourth 
generation with a new CPV phase. 

The NLO PQCD calculation within SM shows @] 
that P' is in the second quadrant, and the color- 
allowed tree amplitude T 1 is roughly real and posi- 
tive. Enhanced by the vertex corrections, C turns 
almost imaginary, and T" + C" is in the fourth quad- 
rant and almost opposite in direction as P' + P^ w ■ As 
Aktt is proportional to the sine of the angle between 
V + C and P' + P EW , its value drops in NLO PQCD 
within SM. In the four generation model at LO in 
PQCD [12j, the amplitude -P EW a-cquires a large CPV 
phase from V t * s Vt'b, which effectively cancels the weak 
phase fa = &rgV* b in T' + C , making Akk° almost 
vanish. The B° — * K + ir~ decay, however, is less af- 
fected, as the -Pew contribution is color-suppressed. 
Since only subleading amplitudes are modified, in 
both proposals the B — > Kir branching ratios are not 
much affected. This is contrary to the proposals based 
on either inelastic [l7| or elastic [HI final state rescat- 
tering, in which leading P' amplitudes are enhanced. 
This is the reason why the B — ► Kir puzzle cannot be 
resolved in such approaches. 

For mixing-induced CP asymmetries, there is 
clearly a major difference between the above pro- 
posals. The time-dependent CP asymmetry of the 
B° — > Ksir mode is defined as 



A CP (B°(t) ^ K s ir°) 
B(B°(t) -» K s ir°) - B{B°(t) 



K s ir°) 



~ B(B°(t) -> A' S 7r ) + B{B°(t) -> K s ir Q ) 

= Ak S tt° cos AmB d t + Sx S ttO sin ArriB d t, (3) 

where Am,B d is the mass difference of the two B- 
meson mass eigenstates, and 



A 



K S TT° 



f 



l + |A Ks7r op 



S 



Ks* 



2 Im Xk s tt° 

1 + \^K s it°\ 2 



(4) 



are the direct and mixing-induced asymmetries, re- 
spectively. The B° — ► Ksir decay has a CP-odd final 
state, and the associated factor, 



A 



a^, M(B° -> Ksir ) 



with M(B° -> K s ir°) = P' - P EW - Ce"*, where 
we have isolated the weak phase in SM [l9[ ■ 



Although C can be enhanced by a few times from 
the vertex corrections in NLO PQCD, within SM 
Skstt is n °t much affected, which stays close to <S ces . 
According to Eq. (O, the leading deviation caused by 
C is proportional to the cosine of the relative strong 
phase between C and P' — P EW - Because the vertex 
corrections also rotate C", it becomes more orthogo- 
nal to P' , and the cosine diminishes. However, in the 
fourth generation model, a NP phase is introduced 
through -P EW , so that A5^- s7r o could be sizable. As we 
will show, the simultaneous accommodation of AAk-k 
and ASx s t T o is nontrivial, with all experimental con- 
straints suitably satisfied. 

The NLO PQCD formalism for the B -> PP, PV 
decays, in which the contributions from the NLO evo- 
lution of Wilson coefficients, the vertex corrections, 
the quark loops, and the magnetic penguin are taken 
into account, can be found in Refs. H and [l6j]. The 
main assumption involved is that we have neglected 
NLO corrections to B meson transition form factors, 
which is not essential for studies of CP asymmetries. 
The procedure to incorporate the fourth generation ef- 
fect can be found in Ref. In Figs.QJa), (b) and (c) 
we plot the predictions for Ak-k, Aktt" and <S^ g7r o vs 
4> sb = &rgV t 1 s V t >b, for three values of r sb = |V^ a Vt'&| = 
0.01, 0.02 and 0.03. We have taken fa = 21.6°, 
fa — 70°, and ra t j = 300 GeV. Theoretical uncer- 
tainties from the above Cabbibo-Kobayashi-Maskawa 
(CKM) matrix elements and from hadronic parame- 
ters, such as Gegenbaucr coefficients in meson distri- 
bution amplitudes, can be included straightforwardly 
as in Ref. [8| . Considering these uncertainties, the pre- 
dicted ranges of A and S will be enlarged by about 
50% and a few percent, respectively. A variation in 
m t i just changes the range of r sb (i.e. the CKM prod- 
uct |V^ s Vt'&|)- Indeed, Ak-k becomes more positive 
for fa, b ~ 90°, while iSr- s7I -o dips downwards, show- 
ing that the nearly vanishing Akk° an d the smaller 
<Sk s k° occur at the same NP phase. Compared with 
Ref. [12] in LO PQCD, we observe that the strength 
of the parameter r sb has shrunk: Ak-k hips sign for 
m t i = 300 GeV and r sb = 0.03 in [12j, but now 
for r sb = 0.02. This implies that the larger color- 
suppressed tree amplitude in NLO PQCD has softened 
the need for larger NP parameters. The decrease of 
R n k, 1.1 for fa. b w 90° shown in Fig. [TJd) is now 
consistent with the observed value. 

The rates and ^Icps of all four B — > Kir modes 
are listed in Table HI which elucidates the effects from 
the NLO Wilson evolution, vertex corrections, quark 
loops, and the magnetic penguin, separately. The ten- 
dency is basically the same as in Ref. @. The main 
change from the incorporation of the fourth generation 
appears in the direct CP asymmetries: Akk° moves 
positive (flips sign) but remains small, and Ak°k° > 
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FIG. 1: (a) Ak-k, (b) A Kv o, (c) S Ks ^o and (d) Ri (i = null, c, n) vs a6 for m t / = 300 GeV with all of the NLO 
corrections. The curves are for r s b = 0.01, 0.02 (not shown in (d)), 0.03, with r s t = 0.03 giving the strongest effect. 



which also depends on -P E w becomes roughly equal 
to Ak-k- Note that Ak-k and Akk° are not in perfect 
match with Eq. (1), but varying hadronic parame- 
ters such as meson distribution amplitudes, the rates, 
Ak-k, and Ak-k could approach data. We do not at- 
tempt any such tuning at present, but took the input 
parameters from earlier NLO PQCD analyses that did 
not involve the fourth generation. 

The predictions for the mixing-induced CP asym- 
metries are relatively insensitive to the hadronic pa- 
rameters. The results for S^Ks and S p ox s are dis- 
played in Fig. O The effect for S^Ks is similar to 
Sk s t T o, going down from S^k s — sin 201 in SM. 
However, S p o Ks behaves in an opposite way: it gets 
enhanced from 0.5 in SM [3, with Eq. © barely 
obeyed. The rates and -4cps of the B — > <pK and 
pK decays are collected in Table [TTJ Compared to the 
NLO results from SM, the branching ratios are not af- 
fected much by the fourth generation as stated before. 
The change mainly appears in the direct CP asymme- 



Quantity +WC +VC +QL +MP +NLO (SM) 



B(A r0 7T ± ) 


23.6 


22.8 


25.2 


17.7 


18.3 (18.1) 


23(^71-°) 


12.6 


12.2 


13.4 


9.6 


9.9 (10.7) 




20.3 


19.8 


21.6 


15.2 


15.8 (15.7) 


B(K°tt°) 


9.3 


9.1 


10.0 


6.8 


7.2 (6.5) 


■Art ■* 





1 








2 (0) 


A KlT o 


-5 


1 


-4 


-7 


1 (-2) 


Ak-k 


-11 


-13 


-9 


-14 


-15 (-12) 


A K o^o 


-5 


-12 


-4 


-7 


-14 (-8) 



TABLE I: Branching ratios (in unit of 10 6 ) and direct 
CP asymmetries (in unit of 10~ 2 ) for B — > Kit decays, 
with m t , = 300 GeV, r sb = 0.025 and 4> ab = 65° [H- The 
labels +WC, +VC, +QL, +MP, and +NLO mean the LO 
results with the NLO Wilson coefficients, the inclusions of 
the vertex corrections, of the quark loops, of the magnetic 
penguin, and of all the above NLO corrections, respec- 
tively. The NLO predictions from SM are presented in the 
parentheses for comparison. 



tries: A^k and A^ro, though small, turn negative; 
ApOfc decreases like Akk°, but is still sizable; A p ox» 
increases like Ak^-k -, and becomes comparable with 
A p r>x and A p k- 

We give a general explanation on the above findings, 
assuming that NP at high energy modifies the Wilson 
coefficients relevant to the electroweak penguin in the 
low-energy effective theory (SM). Both the four gen- 
eration model and the Z' model [2l[ provide such ex- 
amples. The NP effect on the direct CP asymmetries 
and mixing-induced CP asymmetries is written, up to 
corrections of 0(H| W ), as 

AAktt = 2r E wsin(5Ewsin0EW, 
A^A'stt = — 2rEwCOs20i cos^Ewsin^EW, (6) 

where tew is the ratio of the magnitude of the NP 
contribution to the electroweak penguin over the full 
penguin amplitude, <5ew the relative strong phase 
between the electroweak penguin and the full pen- 
guin in SM, and 4>ew the NP phase. Since P E wis 
anti-parallel to X" + C from isospin symmetry [221 ] . 
Jew is expected to be less than 90°. Hence, the NP 
phase 0ew < 90° leads to a decrease of the mag- 
nitude of both Ak-k (if Aktt < in SM, such as 
in the PQCD approach) and Sk s k° . That is, the 
current data of Sk s k° favor the presence of NP in 
the electroweak penguin. Adopting the same parame- 
ters, NP increases S p o Ks as shown in Fig. [2j since the 
NLO PQCD analysis has indicated <5ew > 90° in the 
B — » pK decays [16| , which is attributed to the change 
of the penguin amplitude with destructive combina- 




60 120 180 240 300 360 60 120 180 240 300 360 

0sb <t>sb 

FIG. 2: (a) S^,k s and (b) S p o Ks vs <j> s b with all of the 
NLO corrections. Notation is the same as Fig. [T] 
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Quantity +WC +VC +QL +MP +NLO (SM) 



B((/>K ± ) 


31.0 


12.9 


33.0 


21.7 


8.1 (7.2) 


D\q)A ) 




1 9 n 


OU. 1 


90 9 


i .o yv.o j 


B(p K ) 


r, 1 
0.4 


1.6 


r. o 
o.o 


D.U 


7.7 (7.7) 


Y2{ Z^±\ 

B(p K ) 


o n 
AM 


6. 1 


2.9 


2.9 


o.< (3.9J 




6.0 


7.1 


5.9 


7.0 


7.8 (7.6) 


B(p°K°) 


3.0 


3.7 


2.9 


3.7 


4.2 (4.5) 


A^K 


-3 


-1 


-2 


-5 


-3 (0) 




-3 





-2 


-4 


-2 (2) 


A pK o 


4 


3 


4 


5 


4 (1) 


A p K 


54 


45 


54 


55 


44 (67) 


ApK 


65 


62 


67 


56 


58 (57) 


A p a K o 


19 


27 





16 


27 (6) 



TABLE II: Branching ratios (in unit of 10~ ) and direct 
CP asymmetries (in unit of 10 -2 ) for B — > <f>K and pK 
decays with m t i = 300 GeV, r sb = 0.025 and <f) a b = 65°. 
The notation is the same as Table [I] 



tion of the Wilson coefficients 04 and ciq. If the pat- 



tern of the direct CP asymmetries A p k ~ A p ok 
related to the orientation of the penguin amplitude, 
is confirmed, the increase of S p o Ks will become solid. 
Note that implementing NP into other theoretical ap- 
proaches may not resolve the AAk-k and AiS^^o puz- 
zles, and may not lead to the increase of S p ox s - The 
recent first measurement of S p o Ks [l], Q gives a low 
value, but with rather large errors. 

In this Letter we have studied the effect of a fourth 
generation in the NLO PQCD framework. Combin- 
ing the enhancement of the color-suppressed tree with 
CPV in electroweak penguin renders the predictions 
for AAktt and A5^- s7r o more consistent with data, 
which is a nontrivial result. We predict several other 
mixing dependent and direct CPV effects. Future pre- 
cise measurements at the B factories, the Tevatron, 
the LHC, and (hopefully) Super B factories, will de- 
termine whether NP is called for [23[ by AS < 0, and 
in turn constrain the NP parameters. 
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